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H ! Abstract. Using the general QCD parametrization (GP) we display the magnetic 

moments of the octet baryons including all flavor breaking terms to any order. The 
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hierarchy of the GP parameters allows to estimate a parameter related to the quark 
loops contribution of the proton magnetic moment; its magnitude is predicted to be 
inside a comparatively small interval including the value given recently by Leinweber 



>■ I et al. from a lattice QCD calculation 
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^ ■ The ss contribution /i^ = —^G%j{0) to the proton magnetic moment (compare, for 

^ ; 

the experimental situation, Refs.[Il[2]) has been the object of many QCD calculations 

Oh! 

D . leading to a variety of results. A list of references is found in the table 3 of the review of 

Beck-McKeown [3] and also in recent papers by Leinweber et al. Although, as we 
shall see, our results are compatible with the most recent ones of Leinweber et al. [4], 
it is known that QCD (either in the chiral or in the lattice approaches) has produced 
in the course of the years a variety of predictions for G%j{0), differing in magnitude by 
a factor larger than 10 (and also differing in sign). 

To treat this problem we use here the general QCD parametrization (GP) ([6]- [9]) of 
the magnetic moments of the baryons, including (for the flrst time) terms of all orders 
in flavor breaking. Exploiting a dynamical property of the GP, its "hierarchy" of the 
parameters [ID], we can estimate the order of magnitude of a parameter go related to 
a part of the s quark loops contribution of the proton magnetic moment {go is -30 n 
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in the notation of Ref.[i]). The value of go is predicted inside a rather narrow interval 
that includes the value obtained by Leinweber et al. [1] by a lattice QCD calculationy. 



2. The general QCD parametrization of the magnetic moments 

Consider the operator 

M = (l/2) y'rfV(r X j(r-)) (1) 

where the current 

= eg(x)Q7^g(x) (2) 

is expressed in terms of the quark (m, d, s) fields q{x) and of a charge-diagonal flavor 
matrix Q. The CP parametrizes (compare [6] and [8]-see there in particular footnote 
14) the expectation value of M (Eq.([l])) on the exact eigenstate \iPb) of the QCD 
Hamiltonian for the baryon B at rest as: 

{^Ijb\M\4jb) = {<I>b\V^MV\<I>b) = {WB\Y.9uG,i^,f)\WB) = {Wb\G{ct , f)\WB) (3) 

The symbols in the above equations were used repeatedly in previous work ([6]-dm962), 
but, for convenience, we recall part of them again: Thus \^b) is an auxiliary three 
quarks- no gluon state -having simple properties (e.g. L = 0); is a QCD unitary 
transformation leading from the auxiliary state \^b) to the exact QCD eigenstate \iIjb) 
{\4'b)= V\^b))- In Eq-dS]) the space coordinates are eliminated in the second step, 
as described in |6]; the Gi,'s are operators, linear in the matrix Q, depending only 
on the spin cr and flavor / of three quarks; \Wb) are the usual spin(l/2)-flavor octet 
states (with spin up); the g^'s are parameters independent of Q, the same for the 
magnetic moments {Q = = Diag[2/3, —1/3. — 1/3]) and the Z-moments {Q = 

oc Diag[l - (8/3) sin^ Ow, -1 + (4/3) sin^ Ow, -1 + (4/3) sin^ Ow])- 

^Recently [12], applying the GP to this problem, we related the key formula of Bcck-McKeown [3] 
for /Lip = —{l/3)G%j, to the quark loops corresponding to the Z-trace terms. Note: In our notation 
the negative charge of s is not included in the definition of G%j{0). 

2 



In the limit of exact SU{3) {rriu = md = rris) G(cr, /) is the sum of only three 
possible operators. Writing: Gi = J2iQi(^i, G3 = J^iy^kQi^^k, it is: 

G{(T, f) = + (73G3 + goTr[Q] a, (4) 

i 

In Gi, G3 the index i on Qi arises from the 7 (or Z) being attached to a quark line 



with that index. The meaning o: 



a quark line and its index in the exact state lips) is 



illustrated in the fig.l of Ref.|10j 

The Trace terms (e.g. the last term in Eq.(j4]) and others to appear in the following 
sections) can be due only to quark loops; each loop is connected to the other quark 
lines by gluons; because the quark loops associated to Trace terms including a charge 
must have on their contour an incoming 7 (or Z), the number of gluons connecting 
these loops to the other quark lines must be at least three, due to color conservation 
and the Furry theorem (see, e.g.. Sect. 4 of [H]). From Eq.(jll), the contribution to the 
proton magnetic moment from the s-quarks would be {—1/3)%, but, in the e.m. case, 
this s quark contribution is cancelled by those of the u and d quark loops, due to 

Tr[Q^] = (5) 

The symmetry breaking due to the heavier mass of s (we work here assuming charge 
symmetry = rrid = m) is represented in the GP by terms containing an s projector 

= Diag[0, 0, 1] or the product of two such projectors with different indices (products 
of more than two with different indices do not contribute to matrix elements with 
octet states). Below we list all possible structures. 



3. Flavor breaking in the GP 

i) Terms with one 

■^We recall that in the figure just cited the boxes describe the effect of the transformation V in 
Eq.Q on the three quarks in A quark line may zigzag inside the box as much as one likes, due 

to emission and reabsorption of gluons; also each box contains quark loops and gluon lines connected 
in all possible ways consistent with the prescriptions of QCD. 
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The terms with one to be considered in this subsection were often indicated in 
the past as first order flavor breaking; but it was underlined since the beginning [6]) 
that, because Pf = (-P/)"" for any n, such terms contained also contributions of higher 
order in flavor breaking (Am/m^ = (m^ — m)/ms)- 

There are five operators with one indexed P** (see [HI [8]): They are: G2 = J2i QiPff^i 

; G4 = Y^i^kQiPi^^k ; G5 = Y.ijLkQkPiCTi ; Ge = Y^i^kQiPk^i ; G7 = j^i^^j^^kQiPj'^k- 

Of course they act only on the strange baryons. There is also a Trace term containing 
one P* (that acts, clearly, also on p and n) : 

Go = goTr[QP^]Y,cTi (6) 

i 

The magnetic moments fiB for the octet baryons B are: 

fiB = {^b\{M),\^b) = {Wb\ E gu{G,)z\WB) (7) 

In Eqs. fpi7|) Qq is associated to the flavor breaking s-quark contribution to the proton 
moment, so that 

GUO)=9o + 9o (8) 

(multiplied by —1/3) is the whole s-contribution to fip that one is measuring in the 7-Z 
interference experiments [U [2] . 

At this stage one has 8 parameters and 8 data for the octet (7 magnetic moments 
plus the SA transition moment) but the expectation values of the G^^'s (z/ = 0, . . . 7) 
in the |VrB)'s are not independent (see the identity in Eq.(8) of [9j). Thus yjj it is 
impossible to derive the strangeness content of the nucleon only from the octet-baryons 
magnetic moments without exploiting some additional property of QCD. 

We now list the flavor breaking structures with 2 (or more) P*. 

ii) Terms with two ( or more ) 
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To proceed we consider the following flavor breaking structures with two (they 
act only on S° and 

Gs = J2i^j^k{j>k)QiPjPk^i Gg = Y^i^j-LkQiPkPj^k Gio = J^i^j^k QiPf Pk*^ j 

(9) 

We have also 4 Trace structures of higher order: 

Gg = Tr[QP^] E^ PffT, G[5 = Tr[QP^] E^^k Pf^Tk 

(10) 

Gg = Tr[QP^] E^^k PiP^^k G^ = Tr[QP^] E^^,^kii>,) PfP^fr, 



Using Eq.(l3]) with Q = Q'^, one gets the magnetic moments as linear combinations of 
the parameters Qi^^s. Because the 9 parameters multiplying the structures in Eqs. flQfTOl) 
intervene only in a few combinations Q, it is convenient to define: 

9s = I (-25-76 - 2^70 -98 + 4^9 + fi-io) 

9o=9o + 9S /I IN 

9^ = 9o-9^/^ + H/^ ^ ' 

9^ =90 + 4^?oV3 + 2g'o/3 + 4g'j3 - g^j3 

We get finally (compare Eq. f|T2l) ) the magnetic moments to all orders in flavor breaking 

(in Eq. (fT2l) the baryon symbol indicates the magnetic moment). The contributions of 

the IP* terms (see, e.g. [9J) are obtained by omitting in f[T^ g^, gg and the terms with 

9s, 9o, 9o, 9o defined in the Eq. ffTTl) above: 

P = 91- 9o/3 

n = -{2/3){gi - gs) - go/3 

A = -{l/3){g, -g3 + g2- g^) - g(^/3 

S+ = + (1/9) ((?2 - 4(?4 - 4(75 + 8ge + Sg^) - g^/3 

S- = -{l/3){g, + 2g,) + (l/9)(^2 - ^9a + 2g, - Ag, - 4g,) - g^/3 _ 

2° = -(2/3)((?i - gs) + (2/9)(-2(72 -g^ + 2g, - Ag^ + hg^) +g^- g^/3 

E = -{l/3){gi + 2(73) + (2/9)(-2(72 - g^ - Ag^ - ge - 9? + f^s - Q99) + 9s - g^/3 

/i(SA) = -75(^1 -g3 + 96- 97) 

(12) 

'^On using QiPf = — (1/3)P/ and Tr{QP'^) = —(1/3) many operators in (j9|10|) reduce to the same, 
or to some already considered (e.g. Gq = G2). One may ask why these structures are listed separately; 
the reason is that they correspond to different QCD diagrams; only displaying them separately, the 
order of magnitude of their parameters can be estimated by the "hierarchy"- see Sect. 4. 
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From Eqs. ( fTT][T2l) one gets Eq.( fT3l) . also QCD exact to all orders in flavor breaking. 
Note that /i(SA) can have an imaginary part j6j; in principle this might contribute, via 
other terms, to the S° 7A decay rate. But, because all open S° — A channels are 7 
channels, such terms would imply additional 7's and are therefore negligible. 

V3 /i(SA) = n- (3/2)A - (S+ + S-)/4 + E'-gs + {4/9)g', - {l/9)g^, (13) 



The Eq. f|T3|) (barring the higher order parameters g--, gQ, go) is the Okubo relation |13j . 

From the previous formulas one can calculate the contribution of each specific flavor 

/ = {u, d, s) to the baryon magnetic moments and compare the resulting equations 

with Eqs. (4,5) of Ref.H], obtained exactly from QCD by Leinweber et al. using only 

u, d charge symmetry^ The contribution of the flavor / to the magnetic moment of 

B (except for the loops and to be still multiplied by the charges) is indicated as 

similarly to Ref.jl]. In the GP treatment etc. are calculated as follows: Replace 

the matrix Q in Eq.([2]) with . We then obtain for uP,u", etc. the resultsO 

uP = {A/3)gi + (2/3)^73 
= {-l/3)gi + (4/3)(?3 
= (4/3)(?i + (2/3)^73 + (-2/3)(75 + i^^ge + (4/3)^77 
u^' = i-l/3)gi + (4/3)(73 + (4/3)(75 + h2/3)g, + {A/3)g-r + i-l/3)gs + (4/3)(?9 

(14) 

We easily find using Eqs. fll2|14l) the relations: 

«^^ = S+-S- ; M-° = H°-S- (15) 

uP = 2p + n + go ; = p + 2n + c/o (16) 

that are exact consequences of QCD (plus u, d charge symmetry) since the GP displays 
and takes into account all the operators G^, compatible with QCD. Because of the en- 



■^Charge symmetry for u and d has always been part of our general QCD parametrization, although, 
of course, in some cases -e.g. for the explanation of the extraordinary level of precision of the Coleman- 
Glashow e.m. baryon mass formula [Uj- its applicability was examined in detail (see Sect. 4 of [14]). 

®To obtain (referring to the u quarks) the full contribution /i^ to the magnetic moment of baryon 
B, one must add to the loop contributions -for the proton simply go- and multiply the result by 
the u-charge Qu (for p and n see the Eqs. (20,21) of Ref. ^12j) 
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tirely different derivation, it is expected and satisfactory that Eqs. fll5|16p lead to the 
same equations and therefore to the same constraint (fig. 2 of [4j) between the ratios 
uP jv^ and /u^ as those of [1]. This is seen by inserting in Eqs.f ll6l) the expression 
f fT7|) for qq. To summarize: Because Eqs.f ll5|[T6l) follow exactly from QCD plus charge 
symmetry, the general parametrization developed above -that includes exactly all pos- 
sible terms consistent with QCD, plus charge symmetry- must lead necessarily to the 
above equations. Note, by the way, that some papers (e.g. [15J) misinterpreted the GP 
as some model "a httle more sophisticated than the simplest constituent quark model" . 
The connection of our go^go (Eqs. lHllGllHl) ) with the symbols in [4j is: 

go = -30^. = -G1,(0)(1 (17) 
'i?^ = l + V^o (18) 

4. The hierarchy of the coefRcients in the general QCD parametrization 

Everything, so far, is an exact consequence of QCD (plus u,d charge symmetry). 
Now we turn to the determination of go- As already noted [12] its value cannot be 
determined from the magnetic moments only. One has to exploit a dynamical property 
of QCD, "the hierarchy of the parameters", as anticipated at the end of Sect.l. The 
"hierarchy" is an empirical property, noted long ago [lOj and used in past work (e.g. [161 
[H]); below we recall only its main points. Although the hierarchy was discussed in 
detail in the references cited above, for convenience we will recall briefly (at the end of 
this section) how it works in the first case to which it was applied (the extension of the 
Cell Mann-Okubo mass formula to second order flavor breaking). 

It is a fact that the coefficients g,y of the various structures (the "parameters" ) are 
seen to decrease with increasing complexity of the structure. This was flrst seen in the 
GP of the 8-|-10 baryon masses, but applies generally. It appears from the previous 
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analysis that each baryon property (e.g. the magnetic moments) is QCD parametrized 
as a sum of structures each being a sum of terms having a maximum of three different 
indices (as an example, see the structure Gs = J2i^jytk{j>k)QiPjPk^i)i these indices 
are appended to a quark charge, or to a or to a cr in the term under consideration. 
If the rules of Ref.|10j are adopted to normalize univocally the sums, the coefficients 
of the various structures - that is the parameters Qi, - decrease in a well defined way 
with the number of different indices present in the terms being summed and, of course, 
also with the number of P** (flavor breaking) factors in the structure. Each factor 

n 

P/ produces a reduction [J ^ 1/3 and each pair of indices produces a reduction -as 
discussed in [8]- by a factor from 0.2 to 0.37 (see the footnote 16 of [8J). This "pair 
of indices" reduction is due to the exchange of a gluon between two quark lines and 
applies to baryons or mesons composed of quarks u, d, In the calculations performed 
so far (e.g. [ini [HI [18]) this hierarchy of the coefficients was used to show why certain 
formulas (the Gell Mann-Okubo baryon mass formula, or the Coleman-Glashow or the 
Gal-Scheck formulas) work so well. Observe that to do this -that is to show that the 
terms neglected were indeed negligible- any choice of the "pair of indices" reduction 
factor given above (from 0.2 to 0.37) was equivalent; we selected 0.3 or 0.33 in the 
above evaluations for convenience (the same factor as that due to P^); however, since 
below we are going to use the hierarchy to estimate the order of magnitude of go, we 
will consider the whole interval 0.2 to 0.37 mentioned above. 

To exemplify the above presentation we will summarize the main points of the 
first application of the GP [IE] to the baryon masses. In the general parametrization 
the masses of the eight lowest octet + decuplet baryons can be expressed in terms 
of 8 parameters (the e.m. mass differences are neglected). The known values of the 
^Of course the reduction, for the magnetic moments, is relative to gi « 2.8. 

^Incidentally these reduction factors explain why the naive NRQM, based essentially on terms with 
few indices, works reasonably; to understand this fact on the basis of QCD was precisely the aim of 
the GP (see also [H]). 
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masses allow to calculate all the above parameters; one can check that the parameters 
multiplying the two more complicated structures (sums over 3 indices-containing 2 or 3 
P*) are negligible. Neglecting them one finds a new relation between the masses which 
generalizes the Gell Mann-Okubo formula including to a very good approximation the 
contribution of the 2P* terms. 



5. Estimating go 

For the baryon octet magnetic moments at the IP* level it has been seen repeatedly 
that the hierarchy for the Qi, {v = 1...7) applies fairly well except for the ratio p/n. It is 
well known that for the ratio p/n the hierarchy does not work; indeed the extraordinary 
smallness of [p/n + 3/2) is accidental; therefore one cannot determine the order of 
magnitude of go from the n, n system alone. (For this compare [lO] and the second 
paper in [8]; see also [20])r 

Leaving aside g^, there is however a way to estimate the order of magnitude of go 
from Eqs.f fT2l) : indeed from them one derives the relation (fT9l) (again, of course, an 
exact QCD equation): 

9o = -lp-ln + ^/i(EA) + i(S+ - S-) + 1(S° - S") - + ^9s - ^^?9 (19) 

Inserting experimental values in Eq. (fT9|) one obtains (in nuclear magnetons): 

go = 0.13 ± 0.07 - ^gr + ^gs - ^g9 (20) 

where the error comes essentially frorn |/i(SA)| = 1.61 ± 0.08 (the sign of /i(SA) is 
negative, as can be seen from Eq. (IT^ )J^° 



^Actually the Ml transition A ^ p + 7 allows to understand, via the hierarchy applied to the 
whole octet-decuplet system, the smallness of 173 ; compare the Sects. V, VI -see the Eq.(22)- in [8]. 

^°Neglecting gs and gg, one gets from (|20|) gr = + (2/3)go — 0.09 ± 0.04 instead of gj ~ 0.16, the 
value given in the 2nd paper in [8]. This is because, to determine gr, we now used /i(SA) = — 1.61±0.08 
deduced from its experimental value | /i(SA) |= 1.61 ± 0.08. (Instead in [6], [8] the IP'' calculation 
was used giving ^(SA) = —1.48 ± 0.04/iAr). A new measurement of /^(SA) would be of interest; the 
present value is the result of a Primakoff effect determination [21] dating back to 1986. 
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In ( l20l) and gg are parameters associated to a structure with three indices plus two 
P'^ operators; they are both expected to be of order 0.03. Thus | ((78/12) — {g^/S) \< 
0.01. So the contributions of gs and gg can only affect negligibly the already large 
uncertainty in the first term of Eq. fl20l) . Instead, the parameter g^ is associated to a 
structure with three indices plus one P*; its contribution cannot be neglected. Depend- 
ing on the choice of the gluon exchange reduction factors in the hierarchy (from 0.2 to 
0.37), one expects 

0.036 < I^tI < 0.13 (21) 

We thus obtain from Eq.(l20]), if gj is negative: 0.18 ± 0.07 < < 0.33 ± 0.07 and, if 
c/7 is positive: -0.07 ± 0.07 < < +0.08 ± 0.07. Ahogether: 

- 0.07 ± 0.07 < ^0 < 0.33 ±0.07 (22) 

This interval includes go = 0.28 ± 0.07, the value that can be extracted from Fig.2 and 
Eqs.(ll) of Ref. [4j. As a matter of fact the environment dependence of the u-quark 
magnetic moment as reflected by discussed in |1] and re-derived above using 

the GP with charge symmetry - Eqs. fll5|16l) - suggests to stay in the region of positive 
go, so that the interval to be considered reduces to: 

< c/o < 0.33 ± 0.07 (23) 

Because u"" /u^° = {go + 2n+p)/(E^ — H^) its value is ~ 1.2 for go = 0.33 and ^ 1.7 
for go = 0. As to u^/u^^ = {go + 2p + n)/ (S+ — S^) its values are ~ 1.1 for (70 = 0.33 
and ~ 1.0 for go = 0. 

6. Conclusion 

The main point of this paper has been to show that in G%j{0) =go+go one can 
evaluate go by performing a complete general QCD parametrization of the baryon 
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octet magnetic moments and exploiting the GP "hierarchy of parameters" . The point 
of interest is that by the simple general QCD parametrization one obtains a value 
(compare Eq.( l23i) ) compatible with the value (0.28 ± 0.07) obtained in Ref.[l]. 

Of course ^^^(0) is the sum of Qq and (Eq.([8])), where go is related to 'i?^ of 
Ref.[l] by Eq.(18)p-R^ = 1 + go/go]- Even if 'i?^ had a value two times larger than that 
0.139 ± 0.042 of Ref. [l]-as it seems possible- the value of /if would stay in the region of 
very small values (that is at the limit of the present experimental possibilities) recently 
indicated by HAPPEX (compare the HAPPEX JLAB report 2006 in Ref. [2]) and its 
experimental determination (especially at = 0) would remain extremely difficult. 

[Note added, march 1, 2007] Recently the results of the HAPPEX collaboration 
have been published [22j confirming the above indications. While of course we refer 
to [22] for any detail, we note only that the authors state that the experiments will 
continue at values of higher than the present one 0.1 GeV^. At such low values 
of the results (in the words of the authors) "leave little room for observable nucleon 
strangeness dynamics" . 
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